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ABSTRACT 

Context. The intermediate-mass binary pulsars (IMBPs) are characterized by relatively long spin periods (10 - 200 ms) and massive 
0.4 M Q ) white dwarf (WD) companions. Recently, precise mass measurements have been performed for the pulsar and the WD in 
the IMBP PSR J1802-2124. Some observed properties, such as the low mass of the pulsar, the high mass of the WD, the moderately 
long spin period, and the tight orbit, imply that this system has undergone a peculiar formation mechanism. 
Aims. In this work, we attempt to simulate the detailed evolutionary history of PSR J1802-2124. 

Methods. We propose that a binary system consisting of a neutron star (NS, of mass 1.3 M Q ) and an He star (of mass 1.0 M ), and 
with an initial orbital period of 0.5 d, may have been the progenitor of PSR J1802-2124. Once the He star overflows its Roche lobe, 
He-rich material is transferred onto the NS at a relatively high rate of ~ 10~ 7 - 10~ 6 M Q yr~', which is significantly higher than the 
Eddington accretion rate. A large amount of the transferred material is ejected from the vicinity of the NS by radiation pressure and 
results in the birth of a mildly recycled pulsar. 

Results. Our simulated results are consistent with the observed parameters of PSR J1802-2124. Therefore, we argue that the NS + 
He star evolutionary channel may be responsible for the formation of most IMBPs with orbital periods ;£ 3 d. 
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1. Introduction 

In the standard recycling theory, a neutron star (NS) in a low- 
mass X-ray binary (LMXB) accretes mass from its companion 
star, and the transferred angular momentum causes it to be spun 
up to a short spin period (Alpar et al. 1982; Bhattacharya & van 
den Heuvel 1991). The endpoint of the evolution is a binary con- 
sisting of a millisecond radio pulsar with spin period $10 ms, 
and an He white dwarf (WD). This population is called low- 
mass binary pulsars (LMBPs) (Tauris & van den Heuvel 2006). 
In comparison with LMBPs, another population is that of the 
intermediate-mass binary pulsars (IMBPs), which contain a pul- 
sar with a spin period of tens of milliseconds and a massive CO 
or ONeMg WD with mass of > 0.4 M (Camilo et al. 1996; 
Camilo et al. 2001). 

It is noticeable that the spin periods, magnetic field strengths, 
and orbital eccentricities of IMBPs are considerably greater 
than those of LMBPs (Li 2002), implying that these pulsars are 
mildly recycled. So far, the formation mechanisms of IMBPs 
have not been fully understood. Van den Heuvel (1994) sug- 
gests that these objects originate in low/intermediate-mass X- 
ray binaries (L/IMXBs) with a donor star on the asymptotic 
giant branch. When the donor star overflows its Roche-lobe, 
the unstable mass transfer onto the NS will lead to the for- 
mation of a common envelope (CE). During the spiral-in pro- 
cess, the NS experiences a short timescale, super-Eddington ac- 
cretion and is partially spun up (Ivan den Heuvel & Taam 19841 
Stairs 2004). More recent works indicate that IMBPs are likely 
to evolve through a (sub)thermal timescale mass transfer in 
IMXBs without experiencing CE evolution (King & Ritter 1999; 



Podsiadlowski & Rappaport 2000; Kolb et al. 2000; Tauris et 
al. 2000; Podsiadlowski e t al. 2002; Pfahl et al. 2003. See also 
Pylyser & Savonije 19"88| . |Li (2002)| proposes that the thermal- 



viscous instability of the accretion disks in some IMXBs may 
reduce the amount of mass ttansferred onto the NS and result in 
the birth of a mildly recycled pulsar. 

At present, 16 known IMBPs exist, and only two of them 
have relatively precise measurements of the pulsar masses. For 
PSR J0621 + 1002 the measured mass is M NS = 1-70+JJ j° M 
(INice et al. 20081 The other well-studied IMBP, PSR J1802- 
2124, was discovered by the Parkes Multibeam Pulsar Survey 
(IFaulkner et al. 200"4j l. It has a spin period of 12.6 ms and is in a 
16.8 hr circular orbit with a massive WD companion. By measur- 
ing the general relativistic Shapiro delay, Ferd man et al. (2 010) 
estimate the NS mass to be 1 .24(+0. 11) M and WD mass to be 
0.78(+0.04) M (68% confidence). 

Tauris et al. (2000) have calculated the nonconserva- 
tive evolution of IMXBs and find that IMBPs with a wide 
orbit (P 01 b ^ 3 d) can be formed through a short- 
lived highly super-Eddington mass transfer phase. However, 
this scenario may not work for IMBPs with a short or- 
bital period, such as PSR J 1802-2 12 4, PSRs B0655+64 
( [Jones & Lyne~T988l IHmrner et al. 1995), J1 232-6501, J1435- 
6100 (ICamilo et al. 20011 Manchester 200TI and J1756-5322 
(lEdwards & Bailes 200 H . The CO WD may be formed when 
the donor star evolves to be a He star (van den Heuvel 1994), 
which means the NSs have undergone a short timescale (~ 10 3 
yr) accretion from the donor star in the CE evolution phase 
(Ferdman et al. 2010). However, the detailed evolutionary his- 
tories of IMBPs with a short orbital period like PSR J 1802-2 124 
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have not been explored, so will be investigated in this work. 
Employing the observed parameters of PSR J 1802-2 124, we 
simulate its evolutionary history by using a detailed stellar evo- 
lution code. We argue that a binary system consisting of an NS 
and an He star may be its progenitor. In section 2, we describe 
the formation processes of NS + He star systems and obtain the 
initial parameter space of NS + He star systems. In section 3 
we present a detailed evolutionary path of PSR J 1802-2 124. We 
give a brief summary and discussion in section 4. 



2. Formation of NS + He star systems 

The primordial binaries that formed NS + He star systems con- 
sist of a massive primary and an intermediate-mass secondary 
that produce the NS and the WD, respectively. We used an evo- 
lutionary population synthesis approach based on the rapid bi- 
nary star evolution (BSE) code developed by Hurley et al. (2000, 
2002) to derive the initial parameter space of NS + He star 
systems. For the NS formation, there are two main evolution- 
ary channels: iron core-collapse (ICC) supernovae (SNe) and 
electron-capture (EC) SNe. The former is believed to form a 
higher mass (~ 1.35 M ) population, while the latter produces a 
lower mass (~ 1.25 M ) population (Podsiadlo wski et al. 2005 ; 
Schw ab et al. 20101 . In calculation, we distinguish ICC SNe of 
massive stars from EC SNe or accretion-induced collapse (AIC) 
of ONeMg WDs. Some input parameters including the star for- 
mation rate, the initial mass functions, the CE ejection efficiency 
parameter, the kick velocity dispersion, the distributions of ini- 
tial orbital separations, masses, and the eccentricities are de- 
scribed as followsfl 



1. 



We adopt a constant star formation rate S = 7.6085 yr _1 , 
which implies that one binary with M; > 0.8 M is born in 
the Galaxy per year. 

The primary star mass distribution function is taken to be 
(&(lnMi) = Mi^(Mi), in which the initial mass function 



£(Mi) is given by|Kroupa et al. (1993) 



3. The secondary star mass distribution function is taken to be 
(D(lnM 2 ) = M 2 /Mi = q, which corresponds to a uniform 
distribution of the mass ratio q between and 10. 

4. We take a uniform distribution of lna for the binary separa- 
tion a, namely O(lna) = 0.12328 (Hurley et al. 2002). 

5. The CE ejection efficiency parameter acE, which describes 
the fraction of the released orbital energy used to eject the 
envelope during the CE evolution, i. e. cce = £t>ind/(£orb,f - 
£orb,i) (Hurley et al. 2002). Here £bind i s the total binding 
energy of the primary's envelope, E OI \, f, and £ rb,i are the fi- 
nal and initial orbital energies of the secondary, respectively. 
In this work, we adopt acE = 3 in our standard model and 
&ce = 1 for comparison. 

6. The kick velocity dispersion received by the natal NS. 
For NSs formed by ICC SNe, we adopt relatively stronger 
kick velocities with a dispersion cr ICC = 190 kms~' 
(|Hansen & Phinney 19970 However, for the EC SNe 



1 For other parameters, we refer to Table 3 of Hurley et al. (2002) if 
not mentioned. 

2 Recently, Pinsonneault & Stanek (2006) proposed that two compo- 
nents in 50% of detached binaries have similar masses, i.e. q > 0.87. 

3 Hobbs et al. (2005) find that the observed pulsars could be de- 
scribed by a Maxwellian distribution with a dispersion cr = 265 km s~' . 
However, the population synthesis by Belczynski et al. (2010)| suggests 
that, only when the kick velocity dispersion of natal NS formed in close 
interacting binaries is ~ 170 kms -1 , can the simulated results fit the 



or AIC channel, the explosive energies of electron- 
capture SN explosion are significantly lower than those 
of ICC SNe dDessart et al. 20061 IKitaura et al. 20061 
Jiang et al. 2007). Therefore, it is usually believed 
that the natal NSs formed by both channels may 
receive a smaller kick (IPodsiadlo wski et al. 20041 
IScheck et al. 20041 Poelarends et al. 2008 ; for a review, 
see Ivan den Heuvel 201 01 ). Therefore, we take a moderately 
weak velocity dispersions with ctec.aic = 20 km s _1 Q. 
All stars are assumed to be members of binaries in a circular 
orbit, and they have metallicity of Z = 0.02. In our simula- 
tions, the primary mass Mi, the secondary mass M 2 , and the 
binary separation a range in 0.8 - 80M Q , 0.1 - 80M o , and 
3 - lOOOOfl , respectively. 



First, we show an example of the formation process of an 
NS + He star system. Using the BSE code we calculated the 
evolution of a primordial binary consisting of a primary (star 1) 
of mass Mi = 8.5 M and a secondary (star 2) of mass M 2 
4.4 M (Z = 0.02) in an a = 100.7 R G orbit. At t = 33.11 Myr 
the primary evolves to fill its Roche lobe, and transfers its mate- 
rial to star 2. Because of the relatively high-mass ratio, the mass 
transfer is dynamically unstable, resulting in the formation of a 
CE. After the envelope is ejected, the primary becomes an He 
star (Mi = 1.738 M ), the donor star is still on the main se- 
quence (M 2 = 6.837 M ), but the orbit shrinks to a = 33.19 R Q . 
At t — 40.03 Myr, the He star overflows its Roche lobe once 
more, and the product of the subsequent evolution is a binary 
consisting of an ONeMg WD (Mi = 1 .035 M ) and a main se- 
quence star (M 2 = 7.456 M ), in an orbit of a = 70.915 R e . At 
t = 71.77 Myr, Star 2 starts to overflow its Roche lobe, initi- 
ating the second CE evolution. After the CE is ejected, a close 
binary system consisting of an ONeMg WD (Mi = 1 .035 M ) 
and an He star (M 2 = 1 .459 M Q ) is formed in a tight orbit with 
a = 0.725 R Q (see also |Belczynski & Taam 2004[ ). When the He 
star (the descendant of Star 2) fills its Roche lobe and transfers 
the He-rich material onto the ONeMg WD, the WD may grow to 
the Chandrasekhar limit and collapse to be an NS because of the 
electron capture process (Nomoto 1984) 0. Therefore, through 
two CE phases a compact binary including an NS and a low- 
mass He star (M 2 = 1.049 M ) can be produced. Meanwhile, 
with the sudden increase in the orbital separation caused by the 
received kick during the AIC, the mass transfer ceases. 

To investigate the initial parameter space of NS + He star 
systems, we calculated the evolution of 1 x 10 6 binaries to an age 
of 12 Gyr. Figure 1 shows the distribution of natal NS + He star 
systems in the Me e - ^orb diagram when ctecmc = 20 km s . It 
is clear that most systems have initial orbital periods of P or b,i ~ 
0.01 -Id and initial He star masses M He ,i ~ 0.5 - 3.0 M . For 
the same kick, a lower q-ce results in fewer NS + He star sys- 
tems. This difference develops out of the influence of «ce on the 
evolution of close binaries. A high oqe can prevent coalescence 
during the CE stage, significantly increasing the formation rate 
of NS + He star systems (Liu & Li 2006). 



observed intrinsic ratio between disrupted recycled pulsars and double 
NSs. 

4 For the natal NS via EC SN, |Dewi et al. (2006)1 proposes tnat the 
kick is a Maxwellian distribution with a dispersion of 30 km s -1 . 

5 Here, we present an evolutionary example of an NS + He star sys- 
tem. However, the NS may also be formed directly by an electron- 
capture supernova of the He-star that originated from a close binary 
by Case B mass transfer, and its progenitor star should have a mass in 
the range 8 to 12 M e (Podsiadlowski et al. 2004). 
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Fig. 1. Distribution of initial NS + He star systems in the Mu e - 
P orb diagram when cr ICC = 190 km s" 1 , and cr EC AIC = 20 km s -1 . and the NS mass, respectively. 
The crosses and the open triangles correspond to the CE ejection 
efficiency parameters cce = 3, and 1, respectively. 
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Fig. 2. Evolution of an NS + He star system with Mns.i = 
1 .3 AT©, M He ,i = 1 -0 M Q , and P orbii = 0.5 d. The solid and dashed 
curves denote the evolutionary tracks of the mass transfer rate 
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3. An example evolutionary path to PSR J1 802-21 24 

In this section we show an example of an evolutionary path 
to PSR J1802-2124, to reproduce its observed parameters. 
We used the stellar evolution code developed by Eggleton 
(1971,1972,1973) to calculate the evolutionary sequences of 
an NS + He star system. In the calculations, we adopted 
the stellar OPAE opacities for a low temperature given 
by Rogers & Iglesias (1992) and Alexander & Ferguson (1994) 
and take the ratio of the mixing length to the pressure scale 
height to be 2.0. The overshooting parameter of the He star 
(with a chemical abundance Y = 0.98, Z = 0.02) is set to 
( |Wang et al. 2009) . 

Wlien the He star overflows its Roche lobe and transfers He- 
rich material onto the NS, the maximum accretion rate of the NS 
should be limited by the Eddington accretion rate, which is twice 
that of hydrogen accretion, i. e. M^dd - 3.0 x 10~ 8 M yr _1 . If the 
mass transfer rate -Mne from the He star is greater than M^dd, 
the radiation pressure from the NS will cause the excess mass to 
be lost from the system at a rate M = Mu e + M^dd- In our calcu- 
lations, we consider three types of orbital angular momentum 
loss from the binary system: (1) gravitational-wave radiation, 
(2) magnetic braking: we adopt the induced magnetic braking 
prescription given by Sills et al. (2000), (3) isotropic winds: the 
transferred matter in excess of the Eddington accretion rate is 
assumed to be ejected from the vicinity of the NS in the form 
of isotropic winds and to carry away the specific orbital angular 
momentum of the NS (Tauris & S avonije 1999) . 

Based on the results of Figure 1, we calculated the evolution 
of an NS + He star system with initial masses Mns.i = 1.3 M Q , 
MHe.i = 1 .0 M Q , and initial orbital period P or b,i = 0.5 d. Figures 2 
and 3 show the evolution of the mass transfer rate, the NS mass, 
the orbital period, and the He star mass. At / 17.7 lMyr, the 
He star fills its Roche lobe and transfers mass at a high rate of 
~ 10~ 7 - 10~ 6 M yr" 1 onto the NS, which is significantly higher 
than the Eddington limit. A large amount of the transferred mate- 
rial is ejected by radiation pressure, so the mass of the NS hardly 
gains anything. Because the material is transferred from the less 
massive He star to the more massive NS, the orbital period con- 
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Fig. 3. Evolution of an NS + He star system with Mns.i = 
1 .3 M Q , M He ,i = 1 .0 M Q , and P orbj i = 0.5 d. The solid and dashed 
curves denote the evolutionary tracks of the orbital period and 
the He star mass, respectively. 

tinuously increases. The mass exchange lasts ~ 0.25 Myr, dur- 
ing which the NS has accreted ~ 4% of the transferred mass, 
or ~ 0.01 M Q . This amount of material is, however, sufficient to 
spin the NS up to ^ 16 ms and cause the magnetic field to decay 
(ITaam & van den Heuvel 19861 van den Heuvel 1994). The end- 
point of the evolution is an IMBP consisting a mildly recycled 
pulsar and a CO WD with a mass of 0.8 1 M in an orbital period 
of 0.7 1 d. This evolution provides a plausible path that leads to 
the formation of PSR J 1802-2 124. 



4. Discussion and summary 

At present, the formation channel of IMBPs is not very well 
understood. The recent timing analysis of PSR J 1802-2 124 
provides precise measurements of the pulsar and WD masses 
(Ferdman et al. 2010). The observed properties of this system, 
such as the low mass of the pulsar, the high mass of the WD, 
short orbital period, and moderately long spin period (in com- 
parison with those of LMBPs) imply that this binary may have 
a specific evolutionary history. One possibility is that IMBPs 
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like PSR J1802-2124 may have evolved from IMXBs. Because 
the material is transferred from the more massive donor star to 
the less massive NS in IMXBs, the mass transfer occurs on a 
short thermal timescale. Therefore, Podsia dlowski et al. (2 002) 
suggest that IMXBs spend most of their X-ray active lifetime 
as LMXBs. In this phase some of the binaries would have un- 
dergone unstable accretion owing to an accretion disk instability 
during the LMXB phase, and thus the accretion efficiency was 
significantly decreased (Li 2002). 

In this work, we attempted to explore whether the evolution- 
ary channel of NS + He star system can reproduce the observed 
parameters of PSR J 1802-2 124. We first investigated the initial 
parameter space of NS + He star systems using a population syn- 
thesis approach. Our simulated results show that most NS + He 
star systems have initial orbital periods f or b ; i ~ 0.01 -Id and 
initial He star masses Mu e j ~ 0.5 - 3.0 M . We then performed 
numerical calculations for the evolution of an NS + He star sys- 
tem with M NS ,i = 1.3 M , M He; i = 1-0 M , and P orb ,i = 0.5 d. 
Detailed evolutionary calculations shows that the NS + He star 
evolutionary channel can successfully reproduce the observed 
parameters of PSR J 1802-2 124. The NS + He star systems have 
a compact orbit, therefore we propose that the NS + He star evo- 
lutionary channel may be responsible for the formation of most 
IMBPs with orbital periods $ 3 d. 

Obviously, two main uncertainties exist in our work. The 
first one is the CE evolution. As seen in Section 2, the ini- 
tial parameter space of the NS + He star systems strongly 
depends on the CE ejection efficiency parameters acE, which 
should vary with stellar mass and evolution, but is still poorly 
known. Additionally, some works argue that, an NS may ex- 
perience hypercritical accretion in the CE phase and col- 
lapse to be a black hole (see IChevalier 19931 IBrown 19951 
IBrown et al. 200 11 1. However, it has been found that ~ 60% NS 
with a mass greater than 2.0 M can survive under hypercriti- 



cal accretion ( |Belczynski et al. 2002 1. The second one is the EC 
SNe and AIC process. When the ONeMg core of an asymp- 
totic giant branch star (ISiess 20071 IPoelarends et al. 20 08) or an 
He star (Nom oto 1987b grows to a critical mass, EC SN would 
be triggered. However, the AIC process originates from an 
ONeMg WD accreting from the donor star or from the merger of 
two CO WDs dNomoto & Iben 19851 INomoto & Kondo 19911) . 
Compared with the AIC channel, EC SN can also produce a de- 
tached NS + He star system. To obtain a long enough spin-down 
timescale for the natal NS, our evolutionary scenario favors de- 
tached binary systems. Therefore, it seems that also NS formed 
by EC SN are reasonable progenitors of mildly recycled pulsars. 
Podsiadlowski et al. (2004) suggest that the minimum mass of 
the NS progenitor may be 10 - 12 M for single stars, while this 
value can be 6 - 8 M Q in binaries. The progenitor masses of EC 
SNe are related to the amount of convective overshooting, and 
to the metallicity of the star (Podsiadlowsk i~eTal. 20041 1. A full 
population synthesis for the birthrates of NS + He star systems 
formed by EC SNe and AIC channels is beyond the scope of this 
paper, and will be addressed in future work. 
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